The use of Open Reading frame ESTs (ORESTES) for analysis of the honey bee transcriptome

Abstract Background The ongoing efforts to sequence the honey bee genome require additional initiatives to define its transcriptome. Towards this end, we employed the Open Reading frame ESTs (ORESTES) strategy to generate profiles for the life cycle of Apis mellifera workers. Results Of the 5,021 ORESTES, 35.2% matched with previously deposited Apis ESTs. The analysis of the remaining sequences defined a set of putative orthologs whose majority had their best-match hits with Anopheles and Drosophila genes. CAP3 assembly of the Apis ORESTES with the already existing 15,500 Apis ESTs generated 3,408 contigs. BLASTX comparison of these contigs with protein sets of organisms representing distinct phylogenetic clades revealed a total of 1,629 contigs that Apis mellifera shares with different taxa. Most (41%) represent genes that are in common to all taxa, another 21% are shared between metazoans (Bilateria), and 16% are shared only within the Insecta clade. A set of 23 putative genes presented a best match with human genes, many of which encode factors related to cell signaling/signal transduction. 1,779 contigs (52%) did not match any known sequence. Applying a correction factor deduced from a parallel analysis performed with Drosophila melanogaster ORESTES, we estimate that approximately half of these no-match ESTs contigs (22%) should represent Apis-specific genes. Conclusions The versatile and cost-efficient ORESTES approach produced minilibraries for honey bee life cycle stages. Such information on central gene regions contributes to genome annotation and also lends itself to cross-transcriptome comparisons to reveal evolutionary trends in insect genomes.

Influence of HAART on Alternative Reading Frame Immune Responses over the Course of HIV-1 Infection

Background: Translational errors can result in bypassing of the main viral protein reading frames and the production of alternate reading frame (ARF) or cryptic peptides. Within HIV, there are many such ARFs in both sense and the antisense directions of transcription. These ARFs have the potential to generate immunogenic peptides called cryptic epitopes (CE). Both antiretroviral drug therapy and the immune system exert a mutational pressure on HIV-1. Immune pressure exerted by ARF CD8(+) T cells on the virus has already been observed in vitro. HAART has also been described to select HIV-1 variants for drug escape mutations. Since the mutational pressure exerted on one location of the HIV-1 genome can potentially affect the 3 reading frames, we hypothesized that ARF responses would be affected by this drug pressure in vivo. Methodology/Principal findings: In this study we identified new ARFs derived from sense and antisense transcription of HIV-1. Many of these ARFs are detectable in circulating viral proteins. They are predominantly found in the HIV-1 env nucleotide region. We measured T cell responses to 199 HIV-1 CE encoded within 13 sense and 34 antisense HIV-1 ARFs. We were able to observe that these ARF responses are more frequent and of greater magnitude in chronically infected individuals compared to acutely infected patients, and in patients on HAART, the breadth of ARF responses increased. Conclusions/Significance: These results have implications for vaccine design and unveil the existence of potential new epitopes that could be included as vaccine targets.

A comparative transcriptome analysis reveals expression profiles conserved across three Eimeria spp. of domestic fowl and associated with multiple developmental stages

Coccidiosis of the domestic fowl is a worldwide disease caused by seven species of protozoan parasites of the genus Eimeria. The genome of the model species, Eimeria tenella, presents a complexity of 55-60 MB distributed in 14 chromosomes. Relatively few studies have been undertaken to unravel the complexity of the transcriptome of Eimeria parasites. We report here the generation of more than 45,000 open reading frame expressed sequence tag (ORESTES) cDNA reads of E. tenella, Eimeria maxima and Eimeria acervulina, covering several developmental stages: unsporulated oocysts, sporoblastic oocysts, sporulated oocysts, sporozoites and second generation merozoites. All reads were assembled to constitute gene indices and submitted to a comprehensive functional annotation pipeline. In the case of E. tenella, we also incorporated publicly available ESTs to generate an integrated body of information. Orthology analyses have identified genes conserved across different apicomplexan parasites, as well as genes restricted to the genus Eimeria. Digital expression profiles obtained from ORESTES/EST countings, submitted to clustering analyses, revealed a high conservation pattern across the three Eimeria spp. Distance trees showed that unsporulated and sporoblastic oocysts constitute a distinct clade in all species, with sporulated oocysts forming a more external branch. This latter stage also shows a close relationship with sporozoites, whereas first and second generation merozoites are more closely related to each other than to sporozoites. The profiles were unambiguously associated with the distinct developmental stages and strongly correlated with the order of the stages in the parasite life cycle. Finally, we present The Eimeria Transcript Database (http://www.coccidia.icb.usp.br/eimeriatdb), a website that provides open access to all sequencing data, annotation and comparative analysis. We expect this repository to represent a useful resource to the Eimeria scientific community, helping to define potential candidates for the development of new strategies to control coccidiosis of the domestic fowl. (C) 2011 Australian Society for Parasitology Inc. Published by Elsevier Ltd. All rights reserved.

Conformational stability of recombinant manganese superoxide dismutase from the filamentous fungus Trichoderma reesei

Superoxide dismutases (SODS; EC 1.15.1.1) are part of the antioxidant system of aerobic organisms and are used as a defense against oxidative injury caused by reactive oxygen species (ROS). The cloning and sequencing of the 788-bp genomic DNA from Trichoderma reesei strain QM9414 (anamorph of Hypocrea jecorina) revealed an open reading frame encoding a protein of 212 amino acid residues, with 65-90% similarity to manganese superoxide dismutase from other filamentous fungi. The TrMnSOD was purified and shown to be stable from 20 to 90 degrees C for 1 h at pH from 8 to 11.5, while maintaining its biological activity. (C) 2011 Elsevier B.V. All rights reserved.

Recombinant expression, localization and in vitro inhibition of midgut cysteine peptidase (Sl-CathL) from sugarcane weevil, Sphenophorus levis

A cDNA coding for a digestive cathepsin L, denominated Sl-CathL, was isolated from a cDNA library of Sphenophorus levis larvae, representing the most abundant EST (10.49%) responsible for proteolysis in the midgut. The open reading frame of 972 bp encodes a preproenzyme similar to midgut cathepsin L-like enzymes in other coleopterans. Recombinant Sl-CathL was expressed in Pichia pastoris, with molecular mass of about 42 kDa. The recombinant protein was catalytically activated at low pH and the mature enzyme of 39 kDa displayed thermal instability and maximal activity at 37 degrees C and pH 6.0. Immunocytochemical analysis revealed Sl-CathL production in the midgut epithelium and secretion from vesicles containing the enzyme into the gut lumen, confirming an important role for this enzyme in the digestion of the insect larvae. The expression profile identified by RT-PCR through the biological cycle indicates that Sl-CathL is mainly produced in larval stages, with peak expression in 30-day-old larvae. At this stage, the enzyme is 1250-fold more expressed than in the pupal fase, in which the lowest expression level is detected. This enzyme is also produced in the adult stage, albeit in lesser abundance, assuming the presence of a different array of enzymes in the digestive system of adults. Tissue-specific analysis revealed that Sl-CathL mRNA synthesis occurs fundamentally in the larval midgut, thereby confirming its function as a digestive enzyme, as detected in immunolocalization assays. The catalytic efficiency of the purified recombinant enzyme was calculated using different substrates (Z-Leu-Arg-AMC, Z-Arg-Arg-AMC and Z-Phe-Arg-AMC) and rSl-CathL exhibited hydrolysis preference for Z-Leu-Arg-AMC (k(cat)/K-m = 37.53 mM S-1), which is similar to other insect cathepsin L-like enzymes. rSl-CathL activity inhibition assays were performed using four recombinant sugarcane cystatins. rSl-CathL was strongly inhibited by recombinant cystatin CaneCPI-4 (K-i = 0.196 nM), indicating that this protease is a potential target for pest control. (C) 2011 Elsevier Ltd. All rights reserved.

Phylogenetic relationship of dengue virus type 3 isolated in Brazil and Paraguay and global evolutionary divergence dynamics

Background: Dengue is the most important mosquito-borne viral disease worldwide. Dengue virus comprises four antigenically related viruses named dengue virus type 1 to 4 (DENV1-4). DENV-3 was re-introduced into the Americas in 1994 causing outbreaks in Nicaragua and Panama. DENV-3 was introduced in Brazil in 2000 and then spread to most of the Brazilian States, reaching the neighboring country, Paraguay in 2002. In this study, we have analyzed the phylogenetic relationship of DENV-3 isolated in Brazil and Paraguay with viruses isolated worldwide. We have also analyzed the evolutionary divergence dynamics of DENV-3 viruses. Results: The entire open reading frame (ORF) of thirteen DENV-3 isolated in Brazil (n = 9) and Paraguay (n = 4) were sequenced for phylogenetic analysis. DENV-3 grouped into three main genotypes (I, II and III). Several internal clades were found within each genotype that we called lineage and sub-lineage. Viruses included in this study belong to genotype III and grouped together with viruses isolated in the Americas within the lineage III. The Brazilian viruses were further segregated into two different sub-lineage, A and B, and the Paraguayan into the sub-lineage B. All three genotypes showed internal grouping. The nucleotide divergence was in average 6.7% for genotypes, 2.7% for lineages and 1.5% for sub-lineages. Phylogenetic trees constructed with any of the protein gene sequences showed the same segregation of the DENV-3 in three genotypes. Conclusion: Our results showed that two groups of DENV-3 genotypes III circulated in Brazil during 2002-2009, suggesting different events of introduction of the virus through different regions of the country. In Paraguay, only one group DENV-3 genotype III is circulating that is very closely related to the Brazilian viruses of sub-lineage B. Different degree of grouping can be observed for DENV-3 and each group showed a characteristic evolutionary divergence. Finally, we have observed that any protein gene sequence can be used to identify the virus genotype.

Identification of unannotated exons of low abundance transcripts in Drosophila melanogaster and cloning of a new serine protease gene upregulated upon injury

Abstract Background The sequencing of the D.melanogaster genome revealed an unexpected small number of genes (~ 14,000) indicating that mechanisms acting on generation of transcript diversity must have played a major role in the evolution of complex metazoans. Among the most extensively used mechanisms that accounts for this diversity is alternative splicing. It is estimated that over 40% of Drosophila protein-coding genes contain one or more alternative exons. A recent transcription map of the Drosophila embryogenesis indicates that 30% of the transcribed regions are unannotated, and that 1/3 of this is estimated as missed or alternative exons of previously characterized protein-coding genes. Therefore, the identification of the variety of expressed transcripts depends on experimental data for its final validation and is continuously being performed using different approaches. We applied the Open Reading Frame Expressed Sequence Tags (ORESTES) methodology, which is capable of generating cDNA data from the central portion of rare transcripts, in order to investigate the presence of hitherto unnanotated regions of Drosophila transcriptome. Results Bioinformatic analysis of 1,303 Drosophila ORESTES clusters identified 68 sequences derived from unannotated regions in the current Drosophila genome version (4.3). Of these, a set of 38 was analysed by polyA+ northern blot hybridization, validating 17 (50%) new exons of low abundance transcripts. For one of these ESTs, we obtained the cDNA encompassing the complete coding sequence of a new serine protease, named SP212. The SP212 gene is part of a serine protease gene cluster located in the chromosome region 88A12-B1. This cluster includes the predicted genes CG9631, CG9649 and CG31326, which were previously identified as up-regulated after immune challenges in genomic-scale microarray analysis. In agreement with the proposal that this locus is co-regulated in response to microorganisms infection, we show here that SP212 is also up-regulated upon injury. Conclusion Using the ORESTES methodology we identified 17 novel exons from low abundance Drosophila transcripts, and through a PCR approach the complete CDS of one of these transcripts was defined. Our results show that the computational identification and manual inspection are not sufficient to annotate a genome in the absence of experimentally derived data.

Cellular interference in craniofrontonasal syndrome: males mosaic for mutations in the X-linked EFNB1 gene are more severely affected than true hemizygotes

Craniofrontonasal syndrome (CFNS), an X-linked disorder caused by loss-of-function mutations of EFNB1, exhibits a paradoxical sex reversal in phenotypic severity: females characteristically have frontonasal dysplasia, craniosynostosis and additional minor malformations, but males are usually more mildly affected with hypertelorism as the only feature. X-inactivation is proposed to explain the more severe outcome in heterozygous females, as this leads to functional mosaicism for cells with differing expression of EPHRIN-B1, generating abnormal tissue boundaries-a process that cannot occur in hemizygous males. Apparently challenging this model, males occasionally present with a more severe female-like CFNS phenotype. We hypothesized that such individuals might be mosaic for EFNB1 mutations and investigated this possibility in multiple tissue samples from six sporadically presenting males. Using denaturing high performance liquid chromatography, massively parallel sequencing and multiplex-ligation-dependent probe amplification (MLPA) to increase sensitivity above standard dideoxy sequencing, we identified mosaic mutations of EFNB1 in all cases, comprising three missense changes, two gene deletions and a novel point mutation within the 5' untranslated region (UTR). Quantification by Pyrosequencing and MLPA demonstrated levels of mutant cells between 15 and 69%. The 5' UTR variant mutates the stop codon of a small upstream open reading frame that, using a dual-luciferase reporter construct, was demonstrated to exacerbate interference with translation of the wild-type protein. These results demonstrate a more severe outcome in mosaic than in constitutionally deficient males in an X-linked dominant disorder and provide further support for the cellular interference mechanism, normally related to X-inactivation in females.

Aberrant transcript produced by a splice donor site deletion in the TECTA gene is associated with autosomal dominant deafness in a Brazilian family

We ascertained a Brazilian family with nine individuals affected by autosomal dominant nonsyndromic sensorineural hearing loss. The bilateral hearing loss affected mainly mid-high frequencies, was apparently stable with an early onset. Microsatellites close to the DFNA8/DFNA12 locus, which harbors the TECTA gene, showed significant multipoint lod scores (32) close to marker D11S4107. Sequencing of the exons and exon-intron boundaries of the TECTA gene in one affected subject revealed the deletion c.5383 + 5delGTGA in the 5' end of intron 16, that includes the last two bases of the donor splice site consensus sequence. This mutation segregates with deafness within the family. To date, 33 different TECTA mutations associated with autossomal dominant hearing loss have been described. Among them is the mutation reported herein, first described by Hildebrand et al. (2011) in a UK family. The audioprofiles from the UK and Brazilian families were similar. In order to investigate the transcripts produced by the mutated allele, we performed cDNA analysis of a lymphoblastoid cell line from an affected heterozygote with the c.5383 + 5delGTGA and a noncarrier from the same family. The analysis allowed us to identify an aberrant transcript with skipping of exon 16, without affecting the reading frame. One of the dominant TECTA mutations already described, a synonymous substitution in exon 16 (c.5331 G<A), was also shown to affect splicing resulting in an aberrant transcript lacking exon 16. Despite the difference in the DNA level, both the synonymous substitution in exon 16 (c.5331 G<A) and the mutation described herein affect splicing of exon 16, leading to its skipping. At the protein level they would have the same effect, an in-frame deletion of 37 amino-acids (p.S1758Y/G1759_N1795del) probably leading to an impaired function of the ZP domain. Thus, like the TECTA missense mutations associated with dominant hearing loss, the c5383 + 5delGTGA mutation does not have an inactivating effect on the protein. (C) 2012 Elsevier B.V. All rights reserved.

Characterization of a thermostable extracellular tannase produced under submerged fermentation by Aspergillus ochraceus

Background: Tannases are enzymes that may be used in different industrial sectors as, for example, food and pharmaceutical. They are obtained mainly from microorganisms, as filamentous fungi. However, the diversity of fungi stays poorly explored for tannase production. In this article, Aspergillus ochraceus is presented as a new source of tannase with interesting features for biotechnological applications. Results: Extracellular tannase production was induced when the fungus was cultured in Khanna medium with tannic acid as carbon source. The extracellular tannase was purified 9-fold with 2% recovery and a single band corresponding to 85 kDa was observed in SDS-PAGE. The native apparent molecular mass was estimated as 112 kDa. Optima of temperature and pH were 40 degrees C and 5.0, respectively. The enzyme was fully stable from 40 degrees C to 60 degrees C during 1 hr. The activity was enhanced by Mn2+ (33-39%) and NH4+ (15%). The purified tannase hydrolyzed tannic acid and methyl gallate with Km of 0.76 mM and 0.72 mM, respectively, and Vmax of 0.92 U/mg protein and 0.68 U/mg protein, respectively. The analysis of a partial sequence of the tannase encoding gene showed an open read frame of 567 bp and a sequence of 199 amino acids were predicted. TLC analysis revealed the presence of gallic acid as a tannic acid hydrolysis product. Conclusion: The extracellular tannase produced by A. ochraceus showed distinctive characteristics such as monomeric structure and activation by Mn2+, suggesting a new kind of fungal tannases with biotechnological potential. Further, it was the first time that a partial gene sequence for A. ochraceus tannase was described.

Genome analysis of DNA repair genes in the alpha proteobacterium Caulobacter crescentus

Abstract Background The integrity of DNA molecules is fundamental for maintaining life. The DNA repair proteins protect organisms against genetic damage, by removal of DNA lesions or helping to tolerate them. DNA repair genes are best known from the gamma-proteobacterium Escherichia coli, which is the most understood bacterial model. However, genome sequencing raises questions regarding uniformity and ubiquity of these DNA repair genes and pathways, reinforcing the need for identifying genes and proteins, which may respond to DNA damage in other bacteria. Results In this study, we employed a bioinformatic approach, to analyse and describe the open reading frames potentially related to DNA repair from the genome of the alpha-proteobacterium Caulobacter crescentus. This was performed by comparison with known DNA repair related genes found in public databases. As expected, although C. crescentus and E. coli bacteria belong to separate phylogenetic groups, many of their DNA repair genes are very similar. However, some important DNA repair genes are absent in the C. crescentus genome and other interesting functionally related gene duplications are present, which do not occur in E. coli. These include DNA ligases, exonuclease III (xthA), endonuclease III (nth), O6-methylguanine-DNA methyltransferase (ada gene), photolyase-like genes, and uracil-DNA-glycosylases. On the other hand, the genes imuA and imuB, which are involved in DNA damage induced mutagenesis, have recently been described in C. crescentus, but are absent in E. coli. Particularly interesting are the potential atypical phylogeny of one of the photolyase genes in alpha-proteobacteria, indicating an origin by horizontal transfer, and the duplication of the Ada orthologs, which have diverse structural configurations, including one that is still unique for C. crescentus. Conclusion The absence and the presence of certain genes are discussed and predictions are made considering the particular aspects of the C. crescentus among other known DNA repair pathways. The observed differences enlarge what is known for DNA repair in the Bacterial world, and provide a useful framework for further experimental studies in this organism.

Modularity and evolutionary constraints in a baculovirus gene regulatory network

Abstract Background The structure of regulatory networks remains an open question in our understanding of complex biological systems. Interactions during complete viral life cycles present unique opportunities to understand how host-parasite network take shape and behave. The Anticarsia gemmatalis multiple nucleopolyhedrovirus (AgMNPV) is a large double-stranded DNA virus, whose genome may encode for 152 open reading frames (ORFs). Here we present the analysis of the ordered cascade of the AgMNPV gene expression. Results We observed an earlier onset of the expression than previously reported for other baculoviruses, especially for genes involved in DNA replication. Most ORFs were expressed at higher levels in a more permissive host cell line. Genes with more than one copy in the genome had distinct expression profiles, which could indicate the acquisition of new functionalities. The transcription gene regulatory network (GRN) for 149 ORFs had a modular topology comprising five communities of highly interconnected nodes that separated key genes that are functionally related on different communities, possibly maximizing redundancy and GRN robustness by compartmentalization of important functions. Core conserved functions showed expression synchronicity, distinct GRN features and significantly less genetic diversity, consistent with evolutionary constraints imposed in key elements of biological systems. This reduced genetic diversity also had a positive correlation with the importance of the gene in our estimated GRN, supporting a relationship between phylogenetic data of baculovirus genes and network features inferred from expression data. We also observed that gene arrangement in overlapping transcripts was conserved among related baculoviruses, suggesting a principle of genome organization. Conclusions Albeit with a reduced number of nodes (149), the AgMNPV GRN had a topology and key characteristics similar to those observed in complex cellular organisms, which indicates that modularity may be a general feature of biological gene regulatory networks.

Modeling and Identification of an Open-frame Underwater Vehicle: The Yaw Motion Dynamics

A semi-autonomous unmanned underwater vehicle (UUV), named LAURS, is being developed at the Laboratory of Sensors and Actuators at the University of Sao Paulo. The vehicle has been designed to provide inspection and intervention capabilities in specific missions of deep water oil fields. In this work, a method of modeling and identification of yaw motion dynamic system model of an open-frame underwater vehicle is presented. Using an on-board low cost magnetic compass sensor the method is based on the utilization of an uncoupled 1-DOF (degree of freedom) dynamic system equation and the application of the integral method which is the classical least squares algorithm applied to the integral form of the dynamic system equations. Experimental trials with the actual vehicle have been performed in a test tank and diving pool. During these experiments, thrusters responsible for yaw motion are driven by sinusoidal voltage signal profiles. An assessment of the feasibility of the method reveals that estimated dynamic system models are more reliable when considering slow and small sinusoidal voltage signal profiles, i.e. with larger periods and with relatively small amplitude and offset.

Experimental model identification of open-frame underwater vehicles

Most of the works published on hydrodynamic parameter identification of open-frame underwater vehicles focus their attention almost exclusively on good coherence between simulated and measured responses, giving less importance to the determination of “actual values” for hydrodynamic parameters. To gain insight into hydrodynamic parameter experimental identification of open-frame underwater vehicles, an experimental identification procedure is proposed here to determine parameters of uncoupled and coupled models. The identification procedure includes: (i) a prior estimation of actual values of the forces/torques applied to the vehicle, (ii) identification of drag parameters from constant velocity tests and (iii) identification of inertia and coupling parameters from oscillatory tests; at this stage, the estimated values of drag parameter obtained in item (ii) are used. The procedure proposed here was used to identify the hydrodynamic parameters of LAURS—an unmanned underwater vehicle developed at the University of São Paulo. The thruster–thruster and thruster–hull interactions and the advance velocity of the vehicle are shown to have a strong impact on the efficiency of thrusters appended to open-frame underwater vehicles, especially for high advance velocities. Results of tests with excitation in 1-DOF and 3-DOF are reported and discussed, showing the feasibility of the developed procedure.

Seeing beyond reading: a survey on visual text analytics

We review recent visualization techniques aimed at supporting tasks that require the analysis of text documents, from approaches targeted at visually summarizing the relevant content of a single document to those aimed at assisting exploratory investigation of whole collections of documents.Techniques are organized considering their target input materialeither single texts or collections of textsand their focus, which may be at displaying content, emphasizing relevant relationships, highlighting the temporal evolution of a document or collection, or helping users to handle results from a query posed to a search engine.We describe the approaches adopted by distinct techniques and briefly review the strategies they employ to obtain meaningful text models, discuss how they extract the information required to produce representative visualizations, the tasks they intend to support and the interaction issues involved, and strengths and limitations. Finally, we show a summary of techniques, highlighting their goals and distinguishing characteristics. We also briefly discuss some open problems and research directions in the fields of visual text mining and text analytics.